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ONLY FINITELY MANY TRIBONACCI DIOPHANTINE 
TRIPLES EXIST 

CLEMENS FUCHS*, CHRISTOPH HUTLE*, NURETTIN IRMAK**, 
FLORIAN LUCA***, AND LASZLO SZALAY**** 


Abstract. Diophantine triples taking values in recurrence sequences have 
recently been studied quite a lot. In particular the question was raised 
whether or not there are finitely many Diophantine triples in the Tribonacci 
sequence. We answer this question here in the affirmative. We prove that 
there are only finitely many triples of integers 1 < u < v < w such that 
uv 1, uw + 1, vw + 1 are Tribonacci numbers. The proof depends on the 
Subspace theorem. 


1. Introduction 

The theory of diophantine tuples has a long history and origins in work by 
Diophantus. The basic question is to construct sets of rationals or integers with 
the property that the product of any two of its distinct elements plus 1 is a 
square. Many facts are known on this problem most notably Dujella’s result (cf. 
m that there are at most finitely many sets of five positive integers with this 
property (that the product of any two distinct elements from the set plus 1 is a 
square); further results and the attractive history of the problem can be found 
in [4] . Besides the main streamline several variations of this problem have been 
investigated. 

The problem of finding bounds on the size m for Diophantine m-tuples with 
values in linear recurrences is one such variation. The first general result is 
due to Fuchs, Luca and Szalay [5] and states that if {un}n>i is a binary recur¬ 
rent sequence satisfying certain conditions, then there are at most finitely many 
triples of positive integers a < b < c such that a6 -|- 1, ac -I- 1 and be + 1 are all 
members of {un}n>i- The sequences of Fibonacci and Lucas numbers satisfy 
the conditions of the above theorem, and all Diophantine triples with values 
in the Fibonacci sequence or in the Lucas sequence were computed in m and 
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[12], respectively. Later, in [13], the method was extended to find all Diophan- 
tine triples with values in a certain parametric family of Lucas sequences which 
includes the balancing sequence as a particular case. 

All the previous works dealt with binary recurrences. Concerning linear re¬ 
currences of higher order we mention |14j and |15) in which it is shown that there 
are no Diophantine quadruples with values in the Tribonacci sequence 
or in the shifted Tribonacci sequence {T„ -|- l}„>i. Recall that this sequence 
has the property that Tq = Ti = 0, T 2 = 1 and T „+3 = r „+2 + Tn+i + Tn for 
all n > 0. In |14] . they conjectured that only finitely many Diophantine triples 
exist with values in the Tribonacci sequence, but by their method it was not 
possible to deal with this problem. 

In this paper, we prove the above conjecture. More precisely, we have the 
following theorem. 

Theorem 1. There are only finitely triples of integers 1 < u < v < w such that 
l + uv = Tx, l + uw = Ty, 1 + vw = Tz (1) 

hold for some positive integers x, y, z. 

First we recall some facts concerning the Tribonacci sequence {r„}„>o. The 
characteristic equation 

has roots a, /3, 7 = /3, where 

1 -I- wi -I- a;2 
“ “ 3 ’ 

and 

uji = Y^19 -I- 3-\/^ and W 2 = ^19 — 3-\/^. 

Further, Binet’s formula is 

T„ = aa"-h 5^”cy" for all n > 0, (2) 

where 

(a-/3)(a-7)’ 

(see m)- Numerically, 

1.83 < 

0.73 < 

0.18 < 

0.35 < 

a " ' Y ijj ^ a 


c = 


{/3 - a){f (7-a)(7-/3) 


= b (3) 


a < 1.84, 

|/3| = |7| = a-'/'< 0.74, 
a < 0.19, 

|6| = |c| < 0.36. 


(4) 


n—3 




n —2 


— X — X — 1 = U 


/3 = 


2 — uji — UI 2 + V3i{uji — UJ 2 ) 


Further 


(5) 
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for all n > 2 (see m)- 

In Section [21 we give an upper bound for the greatest common divisor (gcd) 
of Ty — 1 and — 1. This result plays an important role in the proof of our 
Theorem |TJ Section [3] proves Theorem |T] apart from Lemmas |T] and |21 which are 
proved in Section [d] The main tools in the proof are the Subspace theorem (we 
use the version from [ 6 ]) and the theory of S'-unit equations (cf. [5]), which we 
use to prove Lemma [TJ 

2. On the gcd of Ty - 1 and - 1 
We start with the following: 

Proposition 1. If 4:< y < z, then 

gcd(Ty-l,T,-l) 

Proof. We may assume that y > 5, otherwise j/ = 4, T 4 — 1 = 1, and there is 
nothing to prove. Let k G (0,1) be some constant to be determined later. Let 
d = gcd(ry — 1 , T 2 — 1 ). li y < Kz + 2, then 

d = gcd(ry <Ty-l <Ty < ( 6 ) 

From now on, we assume that y > nz -\- 2. Let \ = z — y<{l — k)z. Then 

— 1 = 0 (mod d) 
aa^ + bP"^ + cy^ — 1 = 0 (mod d). 

Multiplying the second congruence above by and subtracting from the result 
the first congruence above gives 

bpy{a^ — P^) + — y'^) — — 1 ) = 0 (mod d). 

The left-hand side above is just oP{Ty — 1) — (Tz — 1)- It is an algebraic integer 
which is not zero because otherwise we get that — l)/(Ty — 1) G Q, 

which is false for A > 0. Thus, 

dr] = bp^{a^ — P^) + — y^) — {a^ — 1) (7) 

holds with some nonzero algebraic integer y € K := Q(a,P). The Galois group 
G = Gal(Lf/Q) is 53 . Taking norms from IF to Q in Q gives 

d® < (f\NK/Q{r])\ = \NK/Q{dri)\ 

= Y[\a{b)a{py{a{a)^ - a{P)^) 

ctGG 

+ cr(c)cr(y)^(cr(a)^ - cr(y)^) - (cr(a)^ - 1)|. 

Now G has two elements fixing a (the identity and the involution (/3, y)). The 
other four elements of G map a into one of /3, y and one of /3, y into a. Also, 
any element of G induces a permutation of the coefficients a, b, c (by ([3])) in the 
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Binet formula ([2]) and |6| = |c| > a. We study the size of the factors in the 
product in the right-hand side of ([8|): 

(i) If <j{a) = a, then, by the absolute value inequality, 

\a{b)a{fi)y{a{a)^ - cr(/3)^) -f cr(c)cr(7)^(cr(a)^ - - 1)| 

<a\2\bmy) + 2\b\W + a^-l 


= a^ 1 


m 

\a\y/^ 


2\b\ 


a 


z/2 


- 1 1 < 1.3a^ < 


Here, we used the fact that y > 5, so z > 6, and the numerics (jH). 

(ii) If a{a) ^ a, it then follows that, by applying again the absolute value 
inequality, 

\a{b)a{l3)y{a{a)^ - cr(/3)^) -f cr(c)CT(7)^(CT(a)^ - cr(7)^) - {cr{a)^ - 1)| 

= \a{a)\Ty - 1) - (T, - 1)1 < T, - 1 + \a\-^^^iTy - 1) 

< 2T^ < (2a“^)a^ < 0.6a^ 

Here, we used ([5]). 

With inequality (|S]), we get 

d® < (1.3a(i-'=)^)2(0.6a^)^ < 

giving 

rf<Q,(i-«/ 3)7 ( 9 ) 

Taking k = 3/4 to balance between (IS|) and dH), we get d < which is what 

we wanted to prove. □ 


3. Proof of Theorem 1 

We shall assume that there are infinitely many such triples (u,v,w) with 
corresponding (x, y, z) and eventually reach a contradiction. Solving for u, v, w 
in terms of x,y, z from o, we get 

l {T,-i){Ty-i) _ / (r.-i)(r.-i)~ _ l {Ty-i){T,-i) 

V T, -1 ’ Y -1 ’ Y 1 

(10) 

Since 1 < u < u < w, we get 5 < x < y < z. Further, the case z = 7 entails 
(a;,y) = (5,6), but then the corresponding (u,v,w) in (fTUl) are not integers. 
Thus, z > 8. Since u > 1, we get, by ®, 

^x+v-i > - 1) > T, - 1 > - 1 > a^-\ 

giving X + y > z. Hence, y > z/2. Further, 

T.-1|(T,-I)(r,-1). 
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Let di = gcd{Tz - l,Ty - 1) and d 2 = gcd{Tz - l,Tx - 1). So, Tz - 1 \ did 2 - 
Hence, by ([5]) and Proposition [U we have 


x-2 


>Tx>Tx-l>d2> 


Tz-l 


> 


^2 — 3 1 

a — i z/4-4 

I ' 


di a32/4 


> a 


giving 

x> z/4-2> z/l2, (11) 

at least for z > 12. Since there can be only finitely many such triples {x, y, z) 
with z < 12, we can assume x > z/12 for all the triples that we consider from 
now on. 

We now use results from Diophantine approximations. Put 


Oi := —1/a, bi := b/a, Ci := c/a, 


and then 


u = (1 + aia-“ + b^P^a-^ + c^^i^ 

■ (1 + aia"^ + 

• (1 + aia-^ + + ci-f^ 



where x = (x, y, z) and for each i > 1 the numbers du,i are some coefficients 
which are the product between some rational number (actually, a 2-unit) and a 
monomial (with nonnegative degrees in each indeterminate) in ai,&i,ci. Here 
we are taking the real positive root of a and a respectively. The expansion is 
obtained as follows: We use the binomial series to get 

(l-faia-“ -f bi/d^a-^ + 

T 

= I] ) (aia”" + + O , 

where O has the usual meaning (for the error term see e.g. [101 Lemma 2]). Simi- 

1 /2 

larly we expand the terms (1 -I- aio”*' -I- bipya~y + and (1 -f oiq:“^ 

-I- Ci 7 ^a“^) up to index T. By multiplying out the expressions 
we indeed get 

n —1 

2=1 


u = 


+ 0{a-^n, 
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where the integer n depends only on T and the terms are ordered in some 
arbitrary way. Since x > z/12 and x < y < z we obtain 


n—1 



where we have put ||x|| = z = maxja;, y, z}. We shall use this bound later on. 

Observe that and Cu^ii'x.) are linear forms in x with integer 

coefficients in such a way that the coefficients of A„^i(x) are always < 0 and the 
coefficients of and Cu,i{^) are always > 0. 

We would like to use the Subspace Theorem (see e.g. [6]; the version we are 
going to use can also be found in Section 3 of [7] , which’s notation - in particular 
the notion of heights - we will also use) in order to show, that the expansion for 
u is “essentially” finite (compare with [3 [HI [3]; however, observe that the results 
from these sources are not directly applicable and that we have to work through 
the estimates instead). 

For this, we again take K := Q(a,/3) and choose for S the set of all places 
(normalized as usual so that the Product Formula holds; see |S1 IZ]); which are 
either infinite or in the set {u G Mk ■ lalt, ^ 1 V \j5\v ^ 1 V ^ 1}. 

Using the fixed n (depending on T) from above, we will now define n + 1 
linearly independent sets of linear forms with indeterminants {U,Yq, ..., F„_i) 
for each v £ S. 

We have to distinguish two cases according to whether x + y — z is even or odd. 
In other words we have either x + y — z = 2k or x + y — z = 2k+ 1, respectively, 
which can be rewritten as k = (x + y — z)/2 resp. k = {x + y — z — l)/2 with 
k £ Z. We put e £ {0,1} to distinguish between these cases, whereby we put 
e = 0iix + y — z = 2k and e = lifa; + ?/ — z = 2fc + l. By going to a still infinite 
subset of the solutions we may assume that e is fixed. 

For the standard infinite place oo on C, we define 


n — 1 




( 13 ) 
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is satisfied for all vectors 

^ Q,(a:+y-2-e)/2 . q,A„,„_ i (x) (x)^C„,n_i (x) ^ ^ 


(14) 


Observe that it is here where we need the case distinction (e = 0 resp. e = 1) in 
order to assure that the vector y is in Since by our choice of linear forms 

det(lo,?;j ■ • ■ 5 ln,v) = 1 for all places in S and 


nnK)r=«<y>"""t 


V^S 2=0 


this inequality reduces to 


n n < '^(y) 

v£Si—0 

We split up the double product on the left-hand side into 


(15) 


u - Va^yo - \/aW ^ du,iyt 


n 


OO 

V^OO 


|u|„ • 

v&S\Mk,, 


n 

2=0 vGS 


Inserting the vector we see that n«es\MKc=o 1^1'“ — since u is some 

integer. Furthermore, the last double product equals 1 due to the Product 
Formula. Furthermore, we estimate 


n 


\uL < 


v€Mk,o 

v^oo 


V T,-l 


< (aa^ + b/3^ + — l)'^(aa^ -I- 6/3^ -|- — 1)'^ 

( 4 . 4 W)'. 


< 


And finally, the first expression is just 


/ 7 (x-\-y—z — e)l2 

' aa ^ 




A„,i(x) oB„,i(x) C„,i(x) 


7 


i>n 


which, by (HU, is smaller than Cia' 


-T||x||/12 


. Thus 


nni'....(y)U<Ci^a-™'“.(4.4W)l 


v£Si—0 


(16) 
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Now we choose T (and the corresponding n = n{T)) large enough, such that 

MCi < 

Then (ITS)) implies 




< a 26 


-^IIj 
26 11 -^ 


(17) 


v^S i—0 

Inserting the vector (flTl) on the right-hand side of (1151) . we find 

n— 1 

'H(y) < C2 ■ n{u) ■ •H(Q,^«,.(x)^B„,dx)^c„,.(x)) 


2 = 1 


n —1 


<C2-(4-4M)- [| 


„C3||x|| 


< a 


C4||x|| 


using suitable constants. For the second inequality, we used that 72(u) is just the 
product riwgMK 00 l'“k (since u is an integer), which allows a similar estimation 
as above, and also that 

-^(Q,(^+y-2-e)/2)n < ^n||x||^ 

7^(c^A„,.(x)^B„,,(x)^C„.,(x)) < ■H(q,)C’.||x||.^(^)C;3||x||.^(^)C^||x||^ 

by using the maximum of those finitely many expressions for j = 1,..., n. There¬ 
fore 

H(y)"^ > (18) 

Now, because of (113 and (1181) it just remains to pick 6, such that 

a-Sllx|l < c,-5C4||x|| 


holds for all of our infinitely many solutions x. This works, if we choose S < 
r/(26C4). 

Thus we can apply the Subspace Theorem, which says, that all solutions 
(a;, y, z) of (1131) lie in finitely many proper linear subspaces. In particular, there 
are still infinitely many solutions in one of these proper subspaces, giving a finite 
set lu and (new) coefficients and eu,i for i G 1^ in K such that 


u = c ^( 2 ^+ y -^-<^)/2 





ftBu 




dx) 


(19) 


holds for still infinitely many {x,y,z); here we eventually have to go to a still 
infinite subset. A similar argument holds for the other two variables v and w 
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and we get formulas 



( 20 ) 


which also holds for appropriate choice of the coefficients for infinitely many 
of the {x,y,z). So, we have arrived at the situation when there exist finite 
sets of indices /«, ly, Iw and triples of linear forms {Au,i, Bu,i,Cu,i) for i G 
for % G ly aiid Byj^i^ for t G ly]^ such that for 

infinitely many {x,y,z) the quantities u,v,w given by (ITOl) and (1^ fulfill (P). 


Lemma 1. Suppose that (Q]) has infinitely many solutions. Then there exists a 
line in R.^ given hy 

x{t)=rit + Si, y{t)=r2t + S2, - 2 (t)=T 3 < + S 3 

with rationals (ri, r 2 , r 3 , si, S 2 , S 3 ) such that infinitely many of the solutions 
ix,y,z) to m are of the form {x{n),y{n), z(n)) for some integer n. 

We shall give the prove of this lemma in Section |4l 

Assume we have Lemma P Then there are infinitely many n such that 
{x(n),y{n), z(n)) satisfies our equations. Let A be the least common multi¬ 
ple of the denominators of ti , r 2 , rs. Infinitely many of our n will be in the same 
residue class modulo A, let us call it r. Then writing n = mA + r, we get 

(x, y, z) = ((riA)m -b (rri -b si), (r 2 A)m -b (rr 2 -b S 2 ), {r^A)m + {rrs + S 3 )). 

Since r^A are integers for i = 1,2,3, and x,y,z are integers, we get that rri + Si 
are integers for i = 1, 2,3. Thus, replacing n by m, Vi by A and Si by rri + Si 
for i = 1, 2, 3, we may assume that and Si are all integers for i = 1, 2,3. Now 

x + y — z — e_ {ri +r 2 — rfijm ^ Si -b S 2 — S 3 — e 
2 ~ 2 ^ 2 ■ 

Since we have infinitely many m, it follows that infinitely many of them will be 
in the same residue class modulo 2. Say this residue class is d G {0,1}, and write 
m = 2i + 5. Then 

x+y-z-e 

- - -= (ri -b ra - r 3 )£ -b 5, 

where S' G Z or Z -b 1/2. So, now 

u(£) = j (21) 
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is an integer for infinitely many £, where here 

X = x(€) = (x(2£ + S), y(2£ + (5), z{2£ + S)). 

From this, we will now derive a contradiction. 

First we observe, that there are only finitely many solutions of (EID with 
u{£) = 0. That can be shown by using the fact, that a simple non-degenerate 
linear recurrence has only finite zero-multiplicity (see [5] for an explicit bound). 
We will apply this statement here for the linear recurrence in f; it only remains 
to check, that no quotient of two distinct roots of the form 
^C'u,i(x(^)) jg qj; unity or, in other words, that 

^ ^ (22) 

has no solutions in n G Z/{0}, nia < 0 and mp,m^ > 0. Assume relation 
holds. Suppose n is even (if not replace (I22|) by its square). Then we apply the 
complex conjugation automorphism, that invaries a and switches /3 and 7 and 
we get 

(a'"“/3™^7™'’)" = 1. (23) 

Multiplying (1^ and (1^ we get 

(a2™“(/37)'"''+'"^)" = 1. 

Since /Sy = —and n is even, we obtain 

^{2mc-mp-m^)n _ ^ 


and 2 ma — mp — rrij < 0, a contradiction. 

So now, we can assume, that u{£) ^ 0 for still infinitely many solutions. Using 
m, we can write 

{T,-l)u^ = {T,-l){Ty-l) (24) 

and insert the finite expansion in i for u into (1241) . Furthermore, we use 
the Binet formula 0 and write Tx,Ty and as power sums in x,y, and z 
respectively. Using the parametrization (x,y,z) = {rim + si,r 2 m + S 2 ,r 3 m + S 3 ) 
with m = 2 £ or m = 2 £ + 1 as above, we have expansions in i on both sides of 
(IM)) . Since there must be infinitely many solutions in £, the largest terms on 
both sides have to grow with the same rate. In order to find the largest terms, 
we have to distinguish some cases: If we assume, that e„ ^ 0 for infinitely many 
of our solutions, then is the largest term in the expansion of u 

and we have 




= aa aa” 


It follows that = aa'^. The case e„ = 0 for infinitely many of our solutions 
is not possible, because the right-hand side of (|Ml) would grow faster than the 
left-hand side so that (l24)l would be true for at most finitely many of our £. In 
all other cases, we had e„ = y/aa^, where e G {0,1}. But this contradicts the 
following lemma. 
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Lemma 2. a/o ^ K and y/aa ^ K. 

We will prove this lemma also in Section |4l The contradiction proves our 
theorem. □ 

4. The proofs of Lemmas [T] and [2] 

We start with Lemma [21 

Proof of Lemma O We used the Magma Computational Algebra System to 
justify the statement. In Q(a) the coefficient 

1 


a = 


can be written as 


while 


(a-/3)(a-7) 

a 

+ 2q! + 3’ 

„2 


a 


aa = 


+ 2a + 3 

The splitting field K = Q(a, /3) is generated by a zero of the symmetric 
polynomial 

p{x) = + 2x^ — 3x^ + 2x^ — X + 1. (25) 

In order to determine the zeros of the above polynomial, it is easy to see, that 

1 

X -I - = T, 

X 

where t is a or ,5 or 7 . Thus, denoting by e a zero of ()25l) . we obtain 

a + \/a^ — 4 


We have the representations 


e = 


a = —— 2 e'^ + 3e^ — e +1, 

-5e^ + - Ue^ + - 13e + 6 

22 ’ 


further 


and 


a = 


aa = 


-2e^ -3e*+ e-^+ e^+ 8e-2 
22 


Finally, by testing a and aa by the SquareRoot function of Magma, we can 
conclude that neither a nor aa is a square in Q(a, /3). □ 

Let us continue with Lemma |TJ 
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Proof of Lemma\^ Under the assumption that ([T]) has infinitely many solutions, 
we have already deduced that there are finite sets P , Iw and coefficients , Ow 
and for i € ly, j € Iw in K such that 



Since 1 +vw = Tz = aa^ + 6/3^ + 07 *, we get 


aa^ + - a^S = 1 , 


(26) 


where S' is a certain finite sum of monomials. 

This is an S-unit equation with certain coefficients and of a certain length 
where S is the multiplicative group generated by {a,/3, — 1} inside the complex 
numbers (note that 7 = Replacing 7 = equation (E51) is 

of the form 

=0, (27) 

i€l 

where I is some finite set, d are nonzero coefficients and Ki{x) and Ltfx) are 
linear forms in x with integer coefficients. It is assumed of course that if j ^ j 
then {Ki{x.), Li{x.)) ^ {Kj{x.), Lj{'x.)). Assume that i ^ j are in the same non¬ 
degenerate component of a sub-equation of (l27t . Then 

belongs to some finite list of numbers. Since a and /3 are multiplicatively inde¬ 
pendent, if follows that (Ki — Kj){x) = 0 and {Li — Lj){x) = 0. Since at least 
one of these two forms is nonzero, it follows that we may assume that T(x) = 0 
for some non-zero form L. Thus, there exist r^, Si, U for i = 1,2, 3 such that 


X = rip + siq + ti, y = r2P + S2q + t2, z = rsp + Ssq + 

Let A be the least common multiple of the denominators of r^, Si for i = 1, 2,3 
and let po, qo be such that for infinitely many pairs {p,q) we have p = po 
(mod A) and q = qo (mod A). Then p = po + AA, q = qo + A/i, and 

X = (ri A)A -f (siA/i) -H (ripo + Si^o + U) 

y = (r2A)A -I- (s2A^) -|- {r2Po + S2go + h) 

z = (r3A)A -f {soAp) + {r^po + ssqo + to). 

Since r^A, s^A, x, y, z are integers, we conclude that ViPo + Siqo + ti are also 

integers for i = 1, 2,3. Thus, we may assume that in fact r^, s^, are all integers. 
A similar argument now shows that our equation is of the form 

=0, 
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where r = (A, /i), J is some finite set of indices and fi are nonzero coefficients for 
i G J. Again we may assume that (Mi(r), iVi(r)) ^ {Mj{r), Nj{r)) for i ^ j in 
J. Applying again the theorem on non-degenerate solutions to S'-unit equations, 
we get that there exists some finite set of numbers F such that for some i ^ j, 
the quantity 

This shows that each such r lies on a finite collection of lines, and since there 
are infinitely many r to start with, we get that infinitely many of our r are on 
the same line. This finishes the proof of the lemma. □ 
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